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We derive variational expressions for the grand potential or action in terms of the many-body 
Green function G which describes the propagation of particles and the renormalized four-point vertex 
r which describes the scattering of two particles in many-body systems. The main ingredient of the 
variational functionals is a term we denote as the H-functional which plays a role analogously to the 
usual ^-functional studied by Baym (G.Baym, Phys.Rev. 127, 1391 (1962)) in connection with the 
conservation laws in many-body systems. We show that any H-derivable theory is also ^-derivable 
and therefore respects the conservation laws. We further set up a computational scheme to obtain 
accurate total energies from our variational functionals without having to solve computationally 
expensive sets of self-consistent equations. The input of the functional is an approximate Green 
function G and an approximate four-point vertex F obtained at a relatively low computational cost. 
The variational property of the functional guarantees that the error in the total energy is only of 
second order in deviations of the input Green function and vertex from the self-consistent ones that 
make the functional stationary. The functionals that we will consider for practical applications cor- 
respond to infinite order summations of ladder and exchange diagrams and are therefore particularly 
suited for applications to highly correlated systems. Their practical evaluation is discussed in detail. 
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I. INTRODUCTION 

Total energy calculations play an important role in 
condensed matter physics and quantum chemistry. For 
solid state physicists they are essential in predicting 
structural changes and bulk moduli in solids. In chem- 
istry molecular bonding curves and potential energy sur- 
faces are essential to understand phenomena like molec- 
ular dissociation and chemical reactions. However, ac- 
curate total energy calculations are notoriously difficult 
and computationally demanding. In quantum chemistry 
there are advanced wavefunction methods like configura- 
tion interaction and coupled cluster theorji to calculate 
energies but they can only be applied to relatively small 
molecules. In solid state physics most total energy calcu- 
lations for crystals or surfaces are based on density func- 
tional theorji^ where the density functionals are mostly 
based on the local density approximation (LDA) and gen- 
eralized gradient approximations (GGA)'^. These func- 
tionals have had great success but there are many cases 
where the functionals fail, in which case there is no clear 
systematic route to improvement. We have therefore re- 
cently advanced a different scheme which involves varia- 
tional energy functionals of the many-body Green func- 
tion and applied it succesfuUy to calculate total energies 
of atoms, molecule3^i^*S4L& and the electron gas^. A vari- 
ety of such functionals can be systematically constructed 
using diagrammatic perturbation theory in which the dif- 
ferent functionals correspond to different levels of pertur- 
bation theory. For these functionals we use input Green 
functions that are relatively easy to obtain at low compu- 
tational cost, for instance from a local density or Hartree- 
Fock calculation. The variational property of the func- 
tional then assures that the errors in the energy are only 
of second order in the difference between our approxi- 



mate Green function and the actual Green function that 
makes the functional stationary. This is the essential 
feature that allows one to obtain accurate total energies 
at a relatively low computational cost*. The remaining 
question is then how to select approximate variational 
functionals that yield good total energies. 

In a diagrammatic expansion in many-body perturba- 
tion theory the building blocks are the Green function 
line G, which describes the propagation of particles and 
holes, and the interaction line v which in electronic sys- 
tems is represented by the Coulomb repulsion between 
the electrons. From this diagrammatic structure one can 
proceed to construct variational functionals in several 
ways. First of all we can renormalize the Green func- 
tion lines. This leads to a functional that has been stud- 
ied by Luttinger and WardiS and leads to a functional 
we will call the ^-functional $[G', w], depending on the 
dressed Green function and the bare two-particle interac- 
tion V. The Luttinger- Ward functional has been applied, 
with great success, to the calculation of total energies of 
the electron gas^'^*, and atoms and molecules "^i^^SiLSii^. 
This type of functionals has also received considerable 
attention for Hubbard lattice type system3ii*i2iiiiliiSii&. 
Apart from renormalization of the Green function lines, 
we can also decide to renormalize the interaction lines by 
replacing the bare interaction by a dynamically screened 
one, usually denoted by W. This leads to the functional 
^'[G, W] that was investigated in a paper by Hedin^'' and 
later more completely by Almbladh et aL^'^* and which 
has been applied with succes to calculations of the total 
energy of the electron gasSii^ and atoms^. This type 
of functionals has also received considerable attention 
in the Dynamical Mean Field Theory (DMFT) commu- 
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The natural place to use this functional is 



in extended systems in which screening of the long range 
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Coulomb interaction is essential. Finally there is also the 
possibility to renormalize the four-point vertices and re- 
place them by a renormalized four-point vertex F. In this 
work we will concentrate on this type of functionals. 

The natural place for variational functionals of the 
Green function G and the four-point vertex F is in sys- 
tems where short range correlations play an important 
role such as in highly correlated systems. Such a type 
of theory was recently discussed in work of Janisiiii^ 
on the Hubbard model in which it was demonstrated 
how to derive the so-called parquet approximation from 
a functional of the Green function and the four-point 
vertex. Furthermore Katsnelson and Lichtenstein?^^ j^^ve 
considered the electronic structure of correlated metals 
in which the building blocks of the theory are an approx- 
imate T-matrix and a bare or noninteracting Green func- 
tion (or a bare Green function in an effective correlated 
medium when using Dynamical Mean Field Theory ^^). 
For describing the structural properties of such materials 
it would therefore be of great importance to be able to 
calculate the total energy from variational energy expres- 
sions in terms of the Green function and the four-point 
vertex where we use an approximate G and F as an in- 
put. The variational property then guarantees that the 
errors in the energy are only of second order in the devi- 
ations of the input Green function and vertex from the 
true quantities that make the functional stationary. 

The construction of energy functionals in terms of G 
and F is most naturally done by the use of the Hugenholtz 
diagram techniqu92^*2Si2L2& which has the bare four-point 
vertex as a diagrammatic building block. This procedure 
has been carried out in the early 1960s by De Domini- 
(,jg29j30 g^j^jj later in more generality by De Dominicis and 
MartinSi*^ and leads to a functional we will call the El- 
functional S[G, F]. In the latter works the derivation 
has been carried out for a very general many-body sys- 
tem with not only one- and two-body interactions but 
also with i-body and |-body interactions that describe 
Bose-condensed and superconducting phases. Unfortu- 
nately this leads to rather involved equations and dis- 
guises the simpler case in which there are only one and 
two-body interactions. For instance, in the general Bose- 
condensed and super-conducting case no particle-particle 
and particle-hole contributions to the four-point vertex 
can be distinguished. The work of De Dominicis and 
Martin was aimed at demonstrating that one could ex- 
press all thermodynamic quantities completely in terms 
of distribution functions rather than at a practical appli- 
cation of the formalism. In their work there is, therefore, 
no discussion of approximate functionals and of ways of 
evaluating them. However, nowadays the functionals can 
be subjected to numerical computation and it is there- 
fore timely to discuss the formalism from this point of 
view and to present computational schemes to evaluate 
the functionals to calculate total energies. This is exactly 
the purpose of this work. 

If we consider the first of the two papers of De Domini- 
cis and MartinM, we see that they use a purely algebraic 



approach to construct their functional which is not capa- 
ble of displaying its full structure. Their second paper- 
uses a purely diagrammatic approach to derive in much 
more detail the structure of the functional but the deriva- 
tion is quite difficult due to numerous intricate topologi- 
cal theorems that need to be discussed in order to avoid 
double counting of the diagrams. However, we found that 
a combination of both methods discussed in these two 
papers leads to a much quicker derivation of the final re- 
sults. Therefore, in this work we derive, in a as simple as 
possible manner, a variational energy or action functional 
for normal systems using a purely algebraic method in 
combination with a diagrammatic analysis. We use, how- 
ever, one generalization of the formalism of DeDominicis 
and Martin: since the Green functions are generated by 
differentiation of our functionals with respect to time- 
nonlocal potentials, the most natural framework to use 
is the Keldysh Green function techniqueSiSli^Si^. We 
therefore consider generally time-dependent systems that 
are initially in thermodynamic equilibrium. This has two 
other advantages. Firstly it allows for an elegant discus- 
sion of conservation laws which, as was shown by BaymSi, 
are closely connected to $-derivability. Such conserva- 
tion laws were earlier discussed for variational energy 
and action functionals within the 4>- and Vf-formalism in 
connection with time-dependent density-functional the- 
orjiS. In particular we will in this paper show that also 
S-derivable theories are conserving. Secondly, the use of 
finite temperature allows for an elegant treatment of the 
boundary conditions on the Green functions. These are, 
for instance, essential in going from the equations of mo- 
tion for the Green function to the Dyson equation which 
will play an important role in our derivations. 

The paper is divided as follows. We first discuss some 
definitions that form the basis of our subsequent analysis. 
We then derive self-consistent equations that relate the 
Green function and the renormalized four-point vertex. 
Then we provide a general construction of the variational 
functional using purely algebraic methods and we sub- 
sequently analyze the structure of the functional using 
diagrammatic methods. We then briefly discuss the con- 
serving properties of the functional. Finally we discuss 
approximate functionals with details for their practical 
evaluation and present our conclusions and outlook. 



II. DEFINING EQUATIONS 

In the following we will consider a many-body sys- 
tem initially in thermodynamic equilibrium. At an ini- 
tial time Iq the system is subjected to a time-dependent 
external field. The Hamiltonian of the system in an time- 
dependent external potential ly(xt) is (in atomic units) 
given by 

H{t) = h{t) + V (1) 
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where in the usual second quantization notation the one- 
and two-body parts of the Hamiltonian are given by 



h{t) 



(2) 



V ^ -j dxdx'w(r,r')Vi^(x)Vi^(x')?/'(x')^A(x). (3) 

Here x — (r, a) is a space-spin coordinate. The two- 
body interaction wiU usually be taken to be a Coulombic 
repulsion, i.e. w(r,r') = l/|r — r'|. The one-body part of 
the Hamiltonian has the explicit form 



/i(x^) = -^^^ + ^(x^) - 



(4) 



We further introduced the chemical potential /i in the 
one-body part of the Hamiltonian of Eq. I^J in anticipa- 
tion of a finite temperature treatment of the system. We 
first consider the expectation value of an operator O for 
the case that the system is initially in an equilibrium 
state before a certain time t^. For t < to the expectation 
value of operator O in the Schrodinger picture is then 
given by (O) = Tr{pO} where p = e~^^° /TTe~f^"° is the 
density matrix and Hq is the time-independent Hamilto- 
nian that describes the system before the perturbation is 
switched on. We further defined /? = l/ksT, with fc^ the 
Boltzmann constant, to be the inverse temperature, and 
the trace involves a summation over a complete set of 
states in the Hilbert space. After we switch on the field 
the expectation value becomes a time-dependent quan- 
tity given by 



(0)(t)-Tr 



{pOHit)} 



(5) 



where Onit) = U{tQ., t)0{t)U{t, to) is the operator in the 
Heisenberg picture. The evolution operator U of the sys- 
tem is defined as the solution to the equations 



idtU{t,t') = H{t)U{t,t') 
idt'U{t,t') = -U{t,t')H{t') 



(6) 
(7) 



with the boundary condition U{t,t) = 1 . The for- 
mal solution of Eq. © can be obtained by integra- 
tion to yield the time-ordered expression U{t,t') = 
Texp(— j J^/cItH^t)) for t > t' with a similar expres- 
sion with anti-chronological time-ordering for t' > t. The 
operator e~^^° can now be regarded as an evolution op- 
erator in imaginary time, i.e. U{to — iP,to) = e~^^°, if 
we define H{t) to be equal to Ho on the contour running 
straight from to to to — i(3 in the complex time plane. We 
can therefore rewrite our expression for the expectation 
value as 

TT{uito~if3,to)U{to,t)dU{t,to)} 
{0){t) = ^ (8) 



Tr 



{u{to~if3,to)} 



FIG. 1; 
plane 



The Keldysh contour drawn in the complex time 



If we read the time arguments of the evolution operators 
in the numerator of this expression from left to right we 
may say that the system evolves from to along the real 
time axis to t after which the operator O acts. Then 
the system evolves back along the real axis from time 
t to to and finally parallel to the imaginary axis from 
to to to — This observation motivates us to define 
the following action functional (compare with the action 
functionals used in Refsi^L^) 

Y = ilnTT{u{to~iPAo)}, (9) 
where we define the evolution operator on the contour as 

U{to - to) = Tc cxp(-i / dtH{t)). (10) 



Here the integral is taken on the contour and Tc denotes 
time-ordering along the contour'^'*'^^ . When we evaluate 
this quantity for the equilibrium system we see that 



iY 



In Tr 



(3n 



(11) 



where fl is the grand potential. Therefore the to- 
tal energy E of the system is obtained from the zero- 
temperature limit 



iY 

lim — = \im n = E - iiN 



(12) 



where N denotes the number of particles in the system. 
Let us now see how this functional can be used as a gen- 
erating functional by making variations with respect to 
parameters in the Hamiltonian. To do this one needs to 
consider changes in the evolution operator U which are 
readily evaluated using Eqs.® and Q. For instance, 
when we make a perturbation SV{t) in the Hamiltonian 
we have using Eq.@ 

idt SU(t, t') ^ 5V{t)U{t, t') + H{t)dU{t, t') (13) 

with a similar differential equation with respect to t' and 
boundary condition 5tj{t, t) = 0. The solution to this 
equation is given by 

6U{t, t') = ^ dTU{t, t)6V(t)U{t, t') (14) 
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from which variations in the action can be calculated. 
For instance, if we choose the perturbation to be a time- 
dependent and spatially nonlocal potential of the form 

6V{t)^ /dxidx2<5u(xi,X2,i)V^^(xi)7/'(x2) (15) 



we obtain the time-dependent one-particle density matrix 
as a functional derivative with respect to Y 



(^t(xi)^(x2))(t) 



6Y 



(5M(xi,X2,i)' 



(16) 



Similarly, when we consider a time-dependent two-body 
potential of the form 

SV{t) = y"d(XiX2X3X4) 5F(xiX2X3X4, i) 

xV?(xi)^t(x2)^(x3)V'(x4) (17) 

we obtain the time-dependent two-particle density ma- 
trix as a derivative 



(^t(xi)^t(x2)^(x3)V'(x4))(t) 

SY 



(5V^(xiX2X3X4,t) 



(18) 



Note that in order to derive Eas. (|16|l and (|18|l we had 
to make variations 5u and SV for time-variables on the 
contour. After the variation is made all observables are, 
of course, evaluated for physical quantities that are the 
same on the upper and lower branch of the contour. In 
the remainder of the paper we will heavily use the ac- 
tion functional as a generating functional for the many- 
body Green functions. To do this we have to generalize 
the time-local potentials u and V to time-nonlocal ones, 
such that the derivatives of Y with respect to these po- 
tentials become time-ordered expectation values that we 
can identify with the one- and two-particle Green func- 
tions G and G2 ■ By a subsequent Legendre transform we 
then can construct a variational functional in terms of G 
and G2- Let us start out by defining the n-body Green 
function as 

G„(l...n,l'...n') = 

{-iriTciM^) . . . ^PH{n)^pU^') . . . V31,(n')]) (19) 

where we introduced the short notation 1 = (xiti) and 
where we defined the expectation value of a Hcisenberg 
operator as 



Tr 



(O) 



[u{to-if3,to)dHit)} 
Tr{f7(io-*/3,io)} 



(20) 



The many-body Green functions satisfy the following hi- 
erarchy equations^Siii which connect the n-body Green 
function to the n -I- 1 and n — 1 body Green function: 



J2 <5(l/)(-l)"--'G„_i(2 ...n,l'...f-l,j' + l...n') 

—i J (ixi;(xi, x)G'„+i(l . . . rt, x^i, xt]^, 1' . . . n'). (21) 

These equations follow directly from the definition of 
the Green functions, the anti-commutation relations of 
the field operators and the equations of motion of the 
evolution operators in Eq.© and {Tj). The Green func- 
tions are defined for time-arguments on the time contour. 
Such contour Green functions were first introduced by 
KeldyshS and are often denoted as Keldysh Green func- 
|-jQjjg34j35j36 play an important role in nonequilibrium 
systems. The one-particle Green function Gi = G obeys 
the boundary condition G(xito, 2) = —G{xito — i(3, as 
is readily derived using the cyclic property of the trace. 
The property G(l, X2io) — — G(l, X2to ~ */?) for the other 
argument is likewise easily verified as well as similar re- 
lations for the n-body Green functions. These boundary 
conditions are sometimes referred to as the Kubo-Martin- 
Schwinger conditions^Si^ and are essential in solving the 
equations of motion for the Green function^. After 
these preliminaries we are now ready to derive the equa- 
tions that connect the one- and two-body Green functions 
which we will use to construct the variational functional 
Y. 



III. RELATION BETWEEN THE FOUR-POINT 
VERTEX AND THE GREEN FUNCTION 

In order to derive a variational energy functional in 
terms of the Green function G and the renormalized four- 
point vertex T we start out by deriving coupled equations 
between these quantities, similar to the familiar Hedin 
equationsii. However, instead of the usual coupled equa- 
tions in terms of the Green function G and the screened 
interaction W we have equations in terms of the Green 
function G and the four-point vertex F. Since our aim 
is to derive equations in terms of the renormalized four- 
point vertex it is advantageous to write our equations in 
terms of the bare four-point vertex first. This is most 
conveniently done within the Hugenholtz diagram tech- 
nique2S*2SiSL2&. We will therefore first rewrite the two- 
particle interaction as a fourpoint function as 

y ^ \ I dxdx'i;(r,r')?/'^(x)V;^(x')V^(x')V^(x) 



i y d(xiX2X3X4)Vb(xiX2X3X4) 
X (Xl ) -0 1 (X2 ) -0 (X3 ) (X4 ) 



(22) 



where we defined 



Vb(xiX2X3X4) w(ri,r2)[(5(x2 -X3)(5(xi -X4) 

-(5(xi - X3)(5(x2 - X4)] (23) 



{idt,-h{\))Gn{l...n,l' ...n') = 



This term is used as a basic entity in the Hugenholtz 



Vo(1234) = 



G2(1234) = 



2 1 



X 



3 4 



FIG. 2: Four-point vertex corresponding to the Hugenholtz 
diagram technique. 



FIG. 3: Definition of the renormalized four-point vertex F 



diagram technique and is displayed pictorially in fig|21 
We now make use of the fact that the Green function 
can be obtained as a derivative of the functional 



iY{u\ = -lnTr{[/[u](to - i/?,io)} 



(24) 



with respect to a nonlocal (in space and time) potential 
m(12), where 



J7[u](to - i/3,io) = Tcexp(-i j dtH{t)~ 
i I dl I d2tp\xi)u{12)i:{x2)) ■ (25) 



Since this expression contains a double time-integral one 
has to define precisely how the time-ordering in this equa- 
tion is defined. This leads to often overlooked subtleties 
the details of which are presented in Appendix 1X1 where 
we further show that 



G(12) 



. SY[u] 



(26) 



By a subsequent differentiation (see Appendix^) we can 
obtain the two-particle Green function as 



G2(1234) = - 



SGjU) 



G(14)G(23). 



(27) 



If the derivatives are taken at u = we obtain the Green 
functions as defined in Eq. H19|) . If the derivative is taken 
at nonzero u then there is no direct relation between the 
Green function and expectation values of time-ordered 
field operators. However, as shown in the AppendixIXlthe 
Green functions in the presence of a nonlocal potential u 
still satisfy a set of hierarchy equations. The first ones 



{idt, - h{l))G{U') ^ 5{U') 

+ J d2(u(12) + S(12))G(2l') 
(-z9ti-/i(l'))G(ll') = <5(ll') 

d2G(12)(u(2l') + S(21')). 



(28) 



(29) 



where we defined the self-energy operator S and its ad- 
joint S by the equations 



(i2E(12)G(2l') 



d(234)y(1234)G2(432l') 



d2G(12)S(2l') = 



J d(234) G2(1234)F(4321') 



Here we defined 



1/(1234) = v{vi,r2)5{h,t2m23)S{14) 
-<5(13)(5(24)]0i234 



(30) 



(31) 



(32) 



where S{ij) = S{ti,tj)6{xi — Xj) and ^1234 = 1 if ti > 
^2 > ^3 > ^4 (on the contour) and zero otherwise. The 
function 6*1234 therefore ensures that the operators have 
the proper time-ordering before the equal time limits, de- 
scribed by the delta functions, are taken. In the next sec- 
tion we will also allow for more general forms of 1^(1234) 
in order to obtain the two-particle Green function as a 
functional derivative with respect to V. In order to de- 
rive a self-consistent set of equations we have to give a 
relation between the two-particle Green function and the 
self-energy. We first note that S = E. This can derived 
by applying to Eq.(|3TJ the operator {idt^ — h{l)) and to 
Eq.(|2n|) the operator {-idt'^ - h{l')). With the use of the 
equations of motion of the one- and two-particle Green 
functions from Ea. (|21|l the result then follows. As a re- 
mark we note that for more general initial conditions E 
is not longer equal to S^. From the equality of S and 
E it follows that the Green function has a unique inverse 
given by the Dyson equation 



G-\12) 



which satisfies 



{idt, " h{l))6{12) - u(12) - E(12) 
Goi(12) - m(12) - S(12) (33) 



J d2G-^{12)G{21) = J rf2G(12)G-i(2l') =(5(11') 

(34) 

For later reference we also defined the inverse Gg^ of the 
noninteracting Green function in Ea. (|33|l . We now define 
the renormalized four-point vertex by the equation 



G2(1234) = G(14)G(23) - G(13)G(24) 
<i(5678)G(15)G(27)r(5786)G(83)G(64) (35) 



This expression is displayed pictorially in Fig|3| The 
four-point vertex F has the interpretation of a renor- 
malized interaction that describes the scattering of two 
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FIG. 4: Graphical display of the equations that relate the 
selfenergy E to the vertex F. 



particles and will play an important role in our energy 
functional later. By differentiating Eg . (|34|l with respect 
to the nonlocal potential u and using Eas. (|27(l . H3()|l and 
we can readily derive (along the lines of referencoii) 
the equations 

E(18) ^ -i j d(23)V^(1238)G(32) 



- -y (i(234567)y(1234)G(36) 
X G(45)r(5678)G(72) 



(36) 



and 



r(1234) = - 



/ 



d(5678) 



JE(14) 
5G{i2) 
(5S(14) 
(5G(65) 



G(67)G(85)r(7238) (37) 



The Eas. (|SH|l and (|^ . which are pictorially displayed 
in Fig01 represent a self-consistent set of equations 
that generate the perturbation series for the self-energy 
E[G, V] in terms of the Green function and the inter- 
action V . For instance, if one starts by taking F = 
in Ea. (|36|) then from Ea. (|37|) one obtains an improved 
four-point vertex F which inserted in Ea. l|36|l leads to an 
improved self-energy. It should be noted that at a corre- 
lated level the simultaneous solution of these equations 
is in general computationally demanding-'*. Moreover 
the self-consistent vertex and sclf-encrgy are not related 
by a Ward-identity"^^'"*^ although the response functions 
calculated from the nonequilibrium Green function in ex- 
ternal fields still satisfy the /-sumrule^^i^. However, as 
explained in the introduction, our aim is to obtain total 
energies. This can be done without solving the com- 
putationally demanding Eqs.ljSU ^^^^ l(T7|) by using ap- 
proximate Green functions and four-point vertices in the 
variational functional that we will construct in the next 
section. 



IV. CONSTRUCTION OF A VARIATIONAL 
FUNCTIONAL 

In this section we will construct a variational energy 
or action functional of the dressed Green function G and 
the renormalized four-point vertex F. The main reason 



for investigating such a functional is to obtain in a simple 
way contributions to the total energy that correspond to 
the infinite summation of ladder- type diagrams. Such 
diagrams correspond to an infinite number of terms in 
the $ or ^'-functional. In the new variables G and F 
we have a corresponding functional S[G, F] . In order to 
derive the S-functional, which we will denote as the De 
Dominicis functional224^2i2ii224i, we start with the action 
functional 

iY[u,V] = -lnTr{[/[u,l/](io -i/3,io)} (38) 

which we will regard as a functional of u and V , where 
we defined 

U[u, V]{to - if3, to) = Tc cxp{-i I dtHo{t) 



i dl d2i/;t(xi)u(12)?/>(x2) 



-^J d(1234)T/(1234) 

X^t(x,)^t(x2)^(x3)^(x4)) 



(39) 

Here 1/(1234) is a general time-dependent two-body in- 
teraction which we require to have the following symme- 
try properties 

F(1234) = -V^(2134) = -F(1243) = V{2U3) (40) 

This will guarantee that the Feynman rules of the Hugen- 
holtz diagram method are satisfied. Eventually, when we 
have derived our equations, we will set V equal to expres- 
sion Vq of (|32|l . To give precise meaning to expression 
Ea. (|39|) we again have to specify how the time-ordering 
is defined when we expand the exponent. This is done in 
Appendix IbI where we show that 



SY 
^Su{21) 
SY 
'(5^(4321) 



G(12) 



--G2(1234) 



(41) 
(42) 



In Appendix iBlit is further demonstrated that these one- 
and two-particle Green functions, even in the presence 
of time nonlocal fields, are related by the first equation 
of motion of the Martin-Schwinger hierarchy. By a Leg- 
endre transform we can now construct a functional of G 
and G2 

F[G,G2] ^ iY[u[G,G2],V[G,G2]] 



d(12)u(21)G(12) 



-^J d(1234)y(4321)G2(1234) 



(43) 



where we now regard u and V as functionals of G and 
G2. This functional satisfies 



6F 



(5G(12) 
6F 
(5G2(1234) 



-F(4321) 



(44) 
(45) 



Therefore the functional 



iY[G,G2] = F[G,G2] 

d(12)u(21)G(12) 



^^J d(1234)F(4321)G2(1234) 



(46) 



(47) 




Sc = 



for fixed u and is a stationary functional of G and G2 , 
i.e. 



SY 



dG{12) 
. SY 

'(5G2(1234) 



= 

= 



(48) 
(49) 




where we will eventually be interested in the case u = Q 
and V = Vq. We can now modify the functional Y[G, G2] 
such that, rather than the two-particle Green function, 
we can use the renormalized four-point vertex F as a 
basic variable. For this purpose we use Ea. H35|l which is 
displayed pictorially in figEl and which gives G2 [G, F] as 
an explicit functional of G and F. We then define the 
functional 



H[G,T] = F[G, G2[G,F]] 



(50) 



FIG. 5: Graphical display of the self-energy terms. The small 
dot denotes the bare vertex V and the big square denotes the 
full four-point vertex P. 



this term. For fixed G we can also calculate the derivative 
with respect to F for which we have 

SH r SF (5G2(5678) 



(5F(1234) 



= J d(5678) 



(5G2(5678) (5F(1234) 



which is a functional of the Green function G and the 
four-point vertex F. Then for fixed F we have 



-t/(4321) 



(56) 



SH 



SF 



<5G(12) SG{12) 



d(3456) 



5F (5G2(3456) where we defined 



= -it(21) - 11(21) - I]c(21) 
where we defined 



(5G2(3456) <5G(12) 

(51) 



E(12) = E^^(12) + Ec'(12) 



(14) = -i / d(23)y(1234)G(32) 



Sc(18) = 



Sc(18) 



J d(234567)y(1234)G(36) 
xG(45)F(5678)G(72) 

<i(234567)F(1234)G(36) 



xG(45)V'(5678)G(72) 



(52) 
(53) 

(54) 
(55) 



V^(1234) = J d(5678)G(15)G(26)y(5678)G(73)G(84) 

(57) 

which is simply a bare vertex dressed with two ingoing 
and two outgoing dressed Green function lines. Using 
the functional H we can now regard the expression iY of 
Ea. (|47|) as a functional of G and F, i.e. 

iY[G,r] = H[G,r] 

d(12)M(21)G(12) 



(58) 



d(1234)y(4321)G2[G,F](1234) 



From Ea. (|36|l we recognize these terms as selfenergy di- 
agrams. They are displayed graphically in Fig. [S] We 
recognize the first term in Ea. (|53|l for F = Vb as the 
Hartree-Fock part of the self-energy. The second part 
Sc of Eq. (|54|l involving the four-point vertex F describes 
the time-nonlocal correlation part of the self-energy. The 
third part on Ea. H55|l is the adjoint correlation part 
of the self-energy. As mentioned earlier we can show 
from the Kubo-Martin-Schwinger boundary conditions 
for a system initially in thermodynamic equilibrium that 
Sc(12) = Sc(12). However, in the following we will keep 
the tilde on the self-energy to keep track of the origin of 



which is a stationary functional of G and F for fixed u and 
V. We have thus achieved our first goal and expressed 
the action iY as functional of G and F. Our next step is 
to specify the functional H in more detail. The variations 
in H are given by the expression 



SH = 



J d(12)(-u(21) - E(21) - Ec(21))(5G(12) 



-J d(1234)T/(4321)(5F(1234) 
d(12)(G-i(21) - Go 1(21) - Sc(21))5G(12) 

(59) 



-- I d(1234)\/(4321)(5F(1234) 
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2 1 2 1 



G2(1234) = 



X 



4 3 4 3 4 

(a) (6) 

2 1 2 1 ; 



3 4 
(e) (/) 



FIG. 6: Expansion of the 2-particle Green function G2 
terms of the full G. The dot denotes the bare vertex V . 



and hence we see that it is convenient to spht up H as 
follows 

H[G,T] = -tr{ln(-G-i)} 

-tr{Goi(G-Go)}-S[G,r] (60) 

This equation defines a new functional S[G, T] which will 
be the central object for the rest of the paper. In Ea. H60() 
we further defined the trace tr (not to be confused with 
the thermodynamic trace Tr) as 



tvA ^ j dlA{l, r 



(61) 



where 1+ denotes that time ti is approached from above 
on the contour. The definition of the S-functional in 
Ea. (|60|l is convenient because then wc have 



5H 



5G{12) 
SH 
(5r(1234) 



SE 



JG(12) 



SE 



(5r(1234) 



(62) 
(63) 



and therefore from Ea. (|59|l we see that the functional 
S[G,r] satisfies the equations 



SE 



5r(1234) 
_SE_ 

SG{12) 



-^t>(4321) 



Sc(21) 



(64) 
(65) 



The functional E is therefore directly related to the cor- 
relation part of the self-energy. To describe the corre- 
lations in the system it is therefore necessary to further 
study the structure of the S-functional, which we will do 
in detail in the next section. 

Note that in Ea. pUl we could also have written 
ln(G~^) rather than ln(— G~^). These terms differ only 
by a (possibly infinite) constant and depend on the defi- 
nition of the branch cut of the logarithm. However, the 
particular definition here reduces properly to the grand 
potential of the noninteracting system when the interac- 
tions are switched off^^. The final De Dominicis func- 
tional (for M = ) is thus given from Ea. (|58|) and H6U|) 



by 

tY[G, r] = -tr {in(-G-i)} - tr {G^^G - Go)} 
-S[G,r] 

-J J (i(1234)Fo(1234)G2[G,r](4321) (66) 

We can check that in the absence of interactions we have 
iY = — trln{— Gq'^} which yields the grand potential of 
the noninteracting system, as we will discuss in more 
detail later. Let us now check the variational property. 
The derivatives of iY with respect to G and T are given 
by 



SY 



(5r(1234) 
. 6Y 
*(5G(12) 



^(f(4321)-Vb(4321)) 



(67) 



G-i(21)-Goi(21)-Ec(21;y) 
+n2l;Vo) + tci21;Vo) (68) 



where we used that 
i S 



4(5G(56) 



d(1234)Vb(1234)G2[G,r](4321) 

= I](65;Vb) + Sc(65;1/o) (69) 



The variational equations that are obtained by putting 
the derivatives H67() and (|68(l equal to zero, are obviously 
solved for the G and F that self-consistently solve the 
equations 



G-i = Go^ 

Vo = v[G,r] 



S[G,F] 



(70) 
(71) 



where E is calculated from Eas. (|53|) and Ea. (|54|l . There- 
fore the functional y[G,F] is stationary whenever the 
Dyson equation is obeyed and whenever the electron- 
electron interaction expanded in G and F is equal to the 
specified interaction Vq. Equation H66|) for the variational 
functional F[G,F] is the first basic result of this work. 
However, before it can be used in actual calculations we 
have, among others, to specify the specific structure of 
the functional S[G, F]. We will show that for several in- 
finite series of diagrammatic terms contributing to this 
functional we can find explicit expressions in terms of G 
and F. To do this wc first have to study the functional 
V[G, F] of Eg. lfTT)) . This is the topic of the next section. 



V. STRUCTURE OF THE H-FUNCTIONAL 

In this section we analyze in more detail the diagra- 
matic structure of the four-point vertex F and the func- 
tional V^[G, F] of Ea. (|71|l which will allow us to obtain 
more explicit expressions for the functional E. These 
quantities can be directly obtained from a diagrammatic 
expansion of the two-particle Green function. If we ex- 
press the diagrams in terms of the fully dressed Green 
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function G we only need to consider diagrams that do not 
contain any self-energy insertions. Since different authors 
use different definitions and drawing conventions for the 
two-particle Green function, it is important to be clear 
about them. We strictly follow the sign, loop rule and 
drawing conventions of referenceii with the small dif- 
ference that we use Hugenholtz diagram8SSi2SiSLS£. For 
clarity our Feynman rules are given in Appendix [CI In 
fig. El we show the first and second order Hugenholtz di- 
agrams in terms of the fully dressed Green function G 
that contribute to the two-particle Green function G2- 
We see that we can write F as a sum of four classes of di- 
agrams. There are three classes of the form [ab, cd) which 
denote diagrams which by removal of two internal Green 
function lines can separate the diagram into two parts, 
one part being connected to the external points ab and 
one part being connected to points cd. The class (12, 34) 
contains diagrams of the particle-particle type, such as 
diagram (d) in figEJ and will be denoted by Fpp. There 
are two classes of particle-hole type, namely (14, 32) and 
(13, 24) which will be denoted by F^^ and F^^. Examples 
of diagrams of these types are diagrams (e) and (/) in 
figini The remaining diagrams which do not fall into one 
of these classes are denoted by Fq (such as diagram (c) 
in figEll. We can therefore write 

F(1234) = Fj,p(1234) + F^,,(1234) 

+F^J1234)+Fo(1234) (72) 

The simplest diagram in class Fq is simply the bare ver- 
tex 1^(1234) (i.e. diagram (c) in figlHI the factor i follows 
from the Feynman rules in Appendix [2)). Since this dia- 
gram is special we separate it off from Fg and define the 
remaining diagrams Fg by the equation 

Fo(1234) = r;,(1234) + iV{l2M) (73) 

Using Eg. 1)72(1 we can then write 

-iV{l2U) = Fpp(1234) + r^^(1234) 

+F^^(1234) + F;,(1234)-F(1234) (74) 

We will now first show how all the terms on the right- 
hand side of this equation can be constructed as a func- 
tional of F. When we have done this we can insert this 
functional into Eq. H64I) and perform the integration with 
respect to F and thereby construct our desired functional 

Let us start with the particle-particle diagrams Fpp. The 
contribution of all diagrams for Fpp can be written as 
sums of blocks of diagrams J connected with two paral- 
lel Green function lines (see figCI ) . Each of these blocks 
J contains diagrams which cannot disconnect points (12) 
and (34) by cutting two Green function lines (such blocks 
are called simple with respect to (12) and (34) in the ter- 
minology of De Dominicis) and therefore each J-block 
does not contain diagrams of the type Fpp. We thus have 

J(1234) = F(1234) - Fpp(1234) (75) 



1 2 



r(1234) = 




4 3 



FIG. 7: Expression of F in terms of J-blocks 

We introduce a convenient matrix notation 

(12|J|34) = J(1234) (76) 
(12|GG|34) = G(13)G(24) (77) 

Within this notation we can, for instance, conveniently 
write G ~ AB instead of 

G(1234) = (12|AB|34) 

= j (i(56)(12|A|56)(56|B|34) 

= j d(56)A(1256)B(5634) (78) 

If we use this notation, then from the Feynman rules 
in Appendix [0 one can readily convince oneself that in 
matrix notation we simply have 

F — J ~\~ Epp 

= J + ]^JGGJ + {^fJGGJGGJ + . . . 

= J+]^JGGV (79) 

where for every pair of Green function lines we have to 
add a factor of \ (see2SiSL22iH) . This follows because for 
any diagram contributing to J, the diagram with outgo- 
ing lines interchanged leads to the same diagram for F 
(for the simple diagram iV in J it follows from Feynman 
rule 5 in Appendix [G)| . Relation Ea. ((79|l allows us to 
express Fpp in terms of F. We have 

J = F(l + ^GGF)-! (80) 

In combination with Eq.jTSJ this then gives 

Fpp = F - F(l + ^GGF)-! (81) 

which expresses Fpp in terms of F. Let us now do the 
same for the particle-hole diagrams. Since 

Ff;,(1234) = -E;i(2134) (82) 

we only need to construct F^^ as a functional of F. For 
the particle-hole diagrams Ff^ we can follow a similar 
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r(1234)=0 + ^HXH; 



fourpoint function a[r](1234) that we want to integrate 
with respect to T to obtain a functional A, i.e. 



M = y" d(1234)a[r](1234)(5r(1234) 



(90) 



IHXHXH: 



In our case we want to do this for a[T] being Tpp, F^^, 
^ph ^0- Because ST has the symmetry property of 
Ea. H40|l this can also be written as 



FIG. 8: Expression of F in terms of /-blocks 



reasoning as for Tpp and we first write F in terms of 
repeated blocks / given by 

/(1234) = F(1234) - Tf^{123A). (83) 

The expression for F in terms of / is displayed in figlHl 
If we use the notation 



(41|J|23) = 7(1234) 
(12|GG|34} ^ G(31)G(24) 



(84) 
(85) 



where in the first term we defined a new matrix / by 
a cyclic permutation of the indices, then (again using 
the Feynman rules of Appendix Q we have in matrix 
notation 



F = J + F 



ph 



I - IGGI 
I-IGGt 



IGGIGGI 



(86) 



where the alternating signs in Ea. (|86|l are related to 
Feynman rule 4 in Appendix [HI As a remark we note 
that from Ea. ()86|l and (|37|l we see that there is a simple 
relation between / and the self-energy: 



7(1234) = - 



(5E(14) 
(5G(32) 



(87) 



One can indeed check, by iterating equations H36(l and 
(|37|l , that the term (5S /5G only yields diagrams that con- 
tribute to 7. We can now express f in terms of f . We 
have 



which gives 



7 = F(l - GGF)- 



f^, = f - f (1 - GGf )- 



(89) 



Before discussing the last set of diagrams Fq let us see 
if we can integrate the functionals Tpp and F^^ that we 
obtained sofar. To do this we first make a general re- 
mark about functional derivatives. We consider a given 



SA = d(1234) [a(1234) - a(2134) 

-f a(3412) - a(1243)](5F(1234) (91) 

Therefore any part of a which is symmetric in the indices 
(12) or (34) (or anti-symmetric with respect to the inter- 
change of pair (12) and pair (34)) will not contribute to 
this variation. Therefore only certain (anti-)symmetric 
parts of a are uniquely determined as functional deriva- 
tives. This does not pose a problem if the functional a 
we want to integrate already has the same symmetry as 
F. This applies for instance to Tpp and Fq. However, the 
function F^^(1234) is not anti-symmetric in the indices 
(12) and (34). However, the combination 

F^,(1234) - F^,,(2134) = F^,,(1234) + Ff,,(1234) (92) 

has this property and therefore 



2 J d(1234) F^^(1234)(5F(1234) 



d(1234) [F^J1234) + Ff*J1234)]JF(1234) (93) 



We can therefore obtain F^^ - 



- Fp^ as a functional deriva- 



tive by formally integrating F^^ and multiplying the re- 
sulting functional by 2. With this in mind we can now 
address the integration of V in the right hand side of 
Ea. (|64|) with respect to F. Using Eq.fTlJl we can write 



iF(1234) = fpp(1234) 
+ f^,,(1234) 
-I- f;,(1234)- 



Kf^,(1234) 
f(1234) 



(94) 



where the expressions with the tilde are defined as in 
Ea. (|57|) . Let us start by integrating Tpp with respect to 
F. Using Ea. H81|) and taking into account the factor 1/4 
in Ea. H64|l we have 



if. 



GGJ^ ppGG 



4 
1 

4 

SL. 



GGF[1- (l + -GGF)-i]GG 



GGTGG 

pp[G, T] 
ST 



-[I - [l + t:GGT)-^]GG 



(95) 
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where we defined 



Lpp[G,T] = -tr {GGTGGr} - -ti {GGT} 



-tr \ ln(l + -GGT) 



(96) 



In this expression the trace tr (not to be confused with 
the thermodynamic trace Tr) for two-particle functions 
is defined as 



tr{A} = J d(12)(12|A|12) 



(97) 



The diagrammatic expansion of the functional Lpp is dis- 
played in the upper part of figlHI Let us now consider the 
particle-hole diagrams. Since 



tT[GGf^^GGSf]=tT[tf„5T] 



(98) 



it is sufficient to integrate GGTpf^GG with respect to T. 
We have using Eq. (|5^ 



4 pf^ 



icGf [1 - (1 - GGf)-^]GG 

^GGf^^GG + i[l - (1 - GGf)-^]GG 

1 5Lph[G,T] 
2 



5T 



where we defined the functional 



(99) 



Lph[G,T] = itr |GGrGGf} + itr {gGf} 



itr{ln(l-GGf)} 



(100) 



The diagrammatic expansion of the functional Lph is dis- 
played in the lower part of figEI Now since 

tr{(f;^,+f^,J<5r} = 2tr{f^;,5r} 

= 2tr{GGf^^GG(5f} (101) 

we obtain 



^(fph + ^ph) 



5Lph[G,V] 
5V 



(102) 



We now collect our results and define 

L[G,T] = Lpp[G,T]+LpH[G,T] 

-itr {GGTGGT} 
8 

This functional L has the property 



5L _l ,^ 
^ ~ 4 



f^.-f) 



(103) 



(104) 



L„\G,T]-- 



1 1 
3 23 



I 1 
' 42" 




1 1 
+ 5? 




2L,„[G.T]= -- 



FIG. 9: Expansion of the functionals Lpp and Lph in dia- 
grams. The four-point vertex F is denoted with a big black 
dot. 



Using this functional we can now split up the functional 
S further as 



s[G,r] = L[G,r] + L'[G,r] 



(105) 



This defines a new functional L' [G, F] . Then from 
Ea. (|64|) we see that if we differentiate both sides of 
Ea. (|105(l with respect to F we obtain 



,5S i~ l~ 



A 
ph 



ph 



F) + — (106) 



We therefore see by comparing to Ea. (f7l| that the func- 
tional L' must satisfy 



If 



8L'{G,V\ 
5V 



(107) 



This functional can not be written out explicitly, but 
since Fg is well-defined diagrammatically the functional 
h' does have a diagrammatic expansion. The first term 
in this expansion is displayed in fig llOl together with its 
functional derivative. Note that the derivative yields 
four diagrams in accordance with Ea. H91l) . We can 
further consider the functional derivative of the func- 
tional S[G, F] with respect to G. According to Ea. H65() 
this yields self-energy diagrams, as is also seen from 
the diagrammatic expansion of S. The G-derivatives of 
Lpp, Lph and L' lead to correlation self-energy diagrams 
Sc,pp[G,F], Sc,p/i[G,r] and TI^{G,V\ in terms of G and 
F that fall into different topological classes. 

We now again collect our results and find from Eq. H66|l 
that the final De Dominicis functional (for u = ) is given 
by 

zr[G,F] = -tr {in(-G-i)} - tr {G^,\G - Go)} 
-L[G,F] - L'[G,F] - ^tr{VbG2[G,F]} (108) 

We finally write the functional in a different form using 
the Dyson equation of Eq. (|33|) 

iY\G, F] = -tr {ln(S - Gq"^)} - tr {EG} 

-L[G,F] - L'[G,F] - ^tr{FoG2[G,F]} (109) 

We can readily check the variationally property of this 
functional. We then find that 



-tr{[(E-Go-i)"i + G]<5E} 
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L'[GS] = 



5U^ 



\GS\ = - 




FIG. 10: The first term in the expansion of the L'-functional 
and its functional derivative with respect to F. For clarity 
in drawing the diagrams for SL' /ST we interchanged the end- 
point labels rather than making the outgoing lines cross. 



-tr - S(Vb) + Sc(F) - Ec(Vb)]^G} 

+ ^tr|(l/-Vo)(5r} =0 (110) 

whenever y[G, F] = Vq for a self-consistent solution of 
the Dyson equation. The variational functional 11U9I) to- 
gether with the variational property H11U|) is the central 
result of this work. In the next sections we will investi- 
gate the practical evaluation of this functional. It is im- 
portant to note that although the functionals in Ea. (|108l) 
and (|109|l are equivalent when evaluated on the fully self- 
consistent G and F obtained from the Dyson equation 
and F[G,F] = Vq this is not true anymore when evalu- 
ated at approximate G and F. In accordance with ref>^ 
the functional forms in Eas. l|lU8|l and (|109|) will be de- 
noted as the Klein-form and Luttinger- Ward-form of the 
functional Y. It was demonstrated in the ^-formalism 
that the Luttinger- Ward form of the functional is more 
stable (has a smaller second derivative) when used for 
the calculation of total energies^. We will therefore in 
the following use the Luttinger- Ward form of the func- 
tional. 



Then from the variational equation 
SE[GX 



6T 



(111) 



we can construct F[G, Vq] as a functional of G and the 
bare interaction Vq (some examples of this procedure are 
given in the next section). With the functional r[G, Vq] 
defined in this way we define the following $ functional 

<f[G,Vo] = -E[G,r[G,Vo]] - ^tr{VoG2[G,r[G,Vo]]} 

(112) 

and the action functional 

iY[G,Vo] = -tr{ln(S-Goi)} -tr{EG} + $[G,yo] 

(113) 

where in this expression self-energy S[G, r[G, Vq]] must 
also be regarded as a functional of G and Vq. From the 
definition of $ it then follows directly that 

(5$ = -tr|Sc<5G| + ^tr|vb<5r| 

+ tr {(S + Ec)(5g} - ^trjvb'Jr} 



= tr{E^G} 

We therefore obtain the result 
(5$ 



JG(12) 



= S(21) 



(114) 



(115) 



We further have that the functional Y[G, Vq] of Ea. ((TT^ 
is stationary when 



= -tr{((S-Goi)-i + G),5S} 



(116) 



i.e. whenever the Dyson equation is obeyed for a <&- 
derivable self-energy. On the basis of the work of Baym^ 
we can therefore conclude that S-derivable theories are 
conserving. 



VII. APPROXIMATIONS USING THE 
H-FUNCTIONAL 



A. Practical use of the variational property 



VI. H-DERIVABLE THEORIES ARE 
CONSERVING 

In this section we will show that any approximate El- 
functional leads to a corresponding ^-functional. Since 
we know from the work of Baym"^^ that any <I>-derivable 
theory is conserving it follows that also S-derivable the- 
ories are conserving, i.e. they respect the macroscopic 
conservation laws, such a momentum, energy and particle 
number conservation and related constraints such as the 
virial theorem^. Consider any approximate S- functional. 



After having discussed the general properties of the 
functional y[G, F] we will discuss its use in the calcu- 
lation of total energies. For a given approximation to 
S[G,F] the stationary point of the functional Y corre- 
sponds to an approximation for the self-energy and the 
four-point vertex obtained from a solution of the Dyson 
equation and of an equation of Bethe-Salpeter type, both 
of which need to be solved to self-consistency. The so- 
lution of these equations for general electronic systems 
is computationally very expensive or impossible. How- 
ever, if we are only interested in total energies, then we 
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can use the variational property of Y and save greatly in 
computational cost as the full self-consistency step can 
then be skipped. To illustrate this we let G and F be 
self-consistent solutions to the variational equations and 
we let G and F be approximations to G and F. Then we 
have that 



Y[G,r] = 



r s^Y 



1 5G6r 



AGAF 



(117) 



where AG — G ~ G and AF = F — F are the deviations 
from the Green function and four-point vertex to the self- 
consistent ones. We see that the error we make in Y is 
only of second order in AG and AF. We may therefore 
obtain rather accurate energies from rather crude inputs. 
These expectations were indeed borne out by our ear- 
lier calculations within the ^-formalism on atoms and 
molecules®. Obviously the actual error we make also de- 
pends on how large the second derivatives of functional Y 
are. For this reason the Klein and Luttinger-Ward forms 
of the functional perform differently. In fact, experience 
within the ^-functional formalism has shown that the 
Luttinger-Ward is more stable than the Klein functional 
with respect to changes of the input Green functiorJ-. 



B. Approximate .^-functionals 



- j d(234567)Vb(1234)G(36)G(45) 
xVb(5671')G(72) 



(121) 



This amounts to a self-consistent solution of the Dyson 
equation with only second order diagrams. A fully self- 
consistent solution of these equations for molecules was 
recently carried out by us^. One of the next simplest 
approximations is obtained by taking L' = Lp^ = which 
yields the functional 

iYpp[G,T] = -tr{ln(S;-G(7^)} -tr{I]G} 

^tr {GGFGGF} 



-Lpp[G, F] 



--tr{^oG2[G,F]} 

The variational equations correspond to 
G-^ = Gq^ -S[G,F] 



where E is calculated from Eas. (|53|l and Ea. H54|) . 
second variational Ea. H124|l corresponds to 

iVo = T{l + ^GGr)'\ 

This equation can be inverted to give 

r = (iVo){l - lGG{iVo)r^ 



(122) 

(123) 
(124) 

The 
(125) 

(126) 



In the following we study some approximate schemes 
using the S-functional in order to illustrate the formal- 
ism discussed in the preceeding sections. We restrict 
ourselves here to the two most simplest examples, the 
self-consistent second order approximation and the self- 
consistent T-matrix approximation. A more advanced 
approximation, also involving the particle-hole diagrams, 
is discussed in the section on the practical evaluation of 
the S-functional. 

The very simplest nontrivial approximation to the we 
can make to the S-functional is to take Lpp = L 
L' — 0. This yields the functional 



ph 



iY2[G,r 
1 



-tr{ln(i; 
tr {GGFGGF} - 



Go-i)}-tr{EG} 
^tr{VoG2[G,r]} (118) 



and expresses the renormalized four-point vertex as a 
sum of particle-particle (direct and exchange) ladder dia- 
grams in terms of the bare potential Vq. The correspond- 
ing self-energy is then readily obtained from Eas. H53l) and 
(|5Hl by inserting the F of Ea. H126(l in Eq.^^. This ap- 
proximation is equivalent to the self-consistent T-matrix 
approximation. It is clear that the set of approximations 
can be made more and more advanced by using more so- 
phisticated approximations for the S-functional. In the 
following sections we will discuss the numerical evalua- 
tion of iY. We will then among other things, consider 
an approximate four-point vertex obtained from the T- 
matrix approximation as an approximate input for the 
evaluation of the energy functional iY at a more sophis- 
ticated level of perturbation theory. 



which we will denote by Y2 since it only involves second 
order diagrams. The variational equations yield 



G-^ ^ G(7^ -S[G,F] 



1 ~ i ~ 
0-4^-4^° 
which simply implies that T — iVo and that 

E[G,yo](ll') = / d(23)Vb(123l')G(32) 



(119) 
(120) 



VIII. PRACTICAL EVALUATION OF THE 
FUNCTIONAL 

A. Evaluation of the traces 

In this section we discuss the how to evaluate the func- 
tional Y[G, F] in actual applications. Our goal is to eval- 
uate Y for an equilibrium system in which case all two- 
time quantities depend on relative time variables on the 
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vertical stretch of the Keldysh contour. In that case it is 
convenient to go over to a Matsubara representation (we 
use the notation of Kadanoff and Baym^) 



A{z) 



(127) 



A{z) 



dtA{t - t')e"^*^'"> (128) 



where the times are imaginary ( t = —ir for < r < /?) 
and where z = iwR / [3 are the Matsubara frequencies 
which run over even or odd integers n depending on 
whether A is a bosonic or fermionic function. In this 
way the equation of motion for the Green function sim- 
ply attains the form 

(z - /l(xi))G(XiX2,2;) = J(xi,X2) 

+ j rfx3S(xiX3,z)G(x3X2,z) (129) 

For the traces of two-point functions we have the expres- 
sion 

ivA = / dtdyLA{l,l+) 
Jo 

= lim V / dxe''^A(x,x,z) (130) 

r,^0+ ^ J 

For the various traces in the functional Y it is further 
convenient to introduce a one-particle basis, such that 
we can write 

A(X1,X2,Z) = ^ Ay(z)(y5i(xi)(^*(x2) (131) 

Then we have, for instance, that 

trAB = lim ^e''My(z)%(z) (132) 

If we choose the orbitals to be eigenfunctions of the one- 
particle Hamiltonian ft,, 



ft(x)(y5i(x) = ei(pi(x) 



(133) 



then the equation of motion of the Green function attains 
the form 



(z - ei)Gij{z) = 6ij + ^Y.ik{z)Gkj{z) 



(134) 



and we see immediately that the noninteracting Green 
function Go is given by 



Go,,(z) = ^ 



(135) 



Consequently the grand potential for the noninteracting 
system is given by fig = iYq / (3 wherei2*22iiS 



trln{-Goi} 



i z 



(136) 



In the zero-temperature limit (3 oo this simply gives 

N 

e. (137) 



lim rjn = > I 

B^oo ^ 

' 2 = 1 



where the sum runs over the N occupied electron orbitals. 
Note that the chemical potential ^ is included in h such 
that ei = ti — fjL where are the eigenvalues of the one- 
body part of the Hamiltonian. 

As a next step we will discuss how to evaluate the 
functional on an approximate Green function G and an 
approximate vertex F. The input Green function will 
in practice not be a fully interacting Green function but 
rather one obtained from a local density approximation 
(LDA) or from a Hartree-Fock approximation. With ap- 
proximate inputs G and F the first term in Ea. H109|| l can 
be written in a computationally convenient form as^ 

-trln|s[G,f] -G(7^} 

= -tr {in(-G-i)} - tr |ln(l - GSc[G, f ]) j(138) 



where we defined 



I]c[G,F] = E[G,F]-I]«^[G] 



(139) 



and the Green function G from the Dyson equation 

G Go + GoS^^[G]G (140) 

The Green function G therefore presents the first iter- 
ation towards the Hartree-Fock Green function starting 
from G. Therefore G = G"^ when we take G = G^^ 
as an input Green function. The term Ec represents the 
correlation part of the self-energy evaluated at an approx- 
imate G and F. The reason for introducing G is that by 
doing this we have in the last term of Ea. (|138|l elimi- 
nated a static part of the self-energy, which makes this 
term well defined without a convergence factor and also 
makes it decay much faster for large frequencies which 
is computationally advantageous as was shown in RefA . 
The first term in Ea. H138|l can be evaluated analytically 
to give 



iYo = -tr{ln(-G-i)} = - Vln(l 



-He. 



) (141) 



where Ei — ei — fi and are the eigenvalues the Hartree- 
Fock equations with a nonlocal self-energy E^^[G]. In 
practice (for instance for LDA input Green functions) 
these eigenvalues are close to the true Hartree-Fock eigen- 
values. Now the functional Y[G, F] can be written as 

iY[G, f] = iYo 

-tr|ln(l -GEc[G,f])| -tr|l][G,f]G} 

-L[G,f] - L'[G,t] - ^tr{yoG2[G,f]} (142) 



15 



The second term can be evaluated by diagonalization of 
GY^c since for a matrix A{z) we have 

tr{ln(l -^)} Um ^e''Mn(l-Ai(z)) (143) 

where Xi{z) are the eigenvalues of A{z). This completes 
one part of the evaluation of the functional Y . 



B. Evaluation of the L' — 0-functional 

Let us now discuss the evaluation of the L[G,r] and 
L'[G,T] functionals. The evaluation of even the lowest 
order term of the L'-functional will be computationally 
very difficult in practice. The first term in the expansion 
of L' is the pentagon of Fig. (|l()|l containing five fourver- 
tices r. Since every four-point vertex depends on four 
space-time coordinates the pentagon is (apart from the 
spin summations) formally an 80-dimensional integral. 
Fortunately, even the approximation L' = represents 
a very sophisticated many-body approximation. We will 
therefore in the following concentrate on this case and 
consider the evaluation of the functional 

iY[G,r] = iFo - tr {ln(l - GSc[G, r])} -tr{SG} 

-L[G,r]-UY{VoG2[G,r]} (144) 

The evaluation of the first terms in this expression has 
been discussed in the preceding subsection and we will 
therefore concentrate on evaluation of L[G,r]. In this 
term the trace is taken over two-particle functions and 
its evaluation will therefore be slightly different from the 
case discussed above. 

As our approximate T we will take the sum of all 
particle-particle and exchange ladders in terms of Vq for 
which we will eventually take the zero- frequency limit. 
This is the approximate T-matrix used in Ref^. This 
approximate T we will denote as F. This approximate 
F will be expressed in terms of our approximate Green 
function which we will denote with G. Then from Eq. H79() 
we have 



then from Ea. H145|l we see that 7 satisfies 



T = iVo + -VoGGT 



(145) 



If we write 



yo(1234) = S{h,t2)S{h,t4)S{t2,t3)Vo{^l^2^3^4) (146) 

with V()(xiX2X3X4) explicitly given in Ea. lf^ . we see 
that we can write 

f(1234) = 5(ii,i2)(5(i3,i4)7(xiX2X3X4;tit3) (147) 
If we further expand 7 in a basis as 

7(^1X2X3X4; tiia) ^^jijkiitits) 

ijkl 

X V*(xi)^*(x2)(pfc(x3)^,(x4) (148) 



+ 9 X! / d,t2Vo,ijpq 



pqrs 



y.Gqr{tl,t2)Gps{ti,t2)-frskl{t2h) (149) 

which in frequency space attains the form 

li]ki{z) = iVo,yfei - 2^ ^ 51 ^o.'JP9 

z\ pqrs 

xGqr(zi)Gps{z - Zi)jrskl{z) (150) 



(note that for 7 we have to sum over the even Matsubara 
frequencies). For simple approximate Green functions G 
of Hartree-Fock or local density type the frequency sum 
over Zi is readily evaluated. We are now ready to evaluate 
the functionals Lpp[G,T] and Lph[G,T]. They are given 
by the expressions 



Lpp = tr{ln(l + yl)}-tr{A} + -tr{A2} (151) 
Lp,, = itr{ln(l-i?)} + itr{i?} + itr{S2}(152) 



where A = GGF and B = GGF. Therefore in order to 
calculate Lpp and Lph we have to diagonalize A and B 
in a two-particle basis. Let us start by the calculation of 
A. We have for our approximate F and G : 



Aijki{tit2hti) = (^(^3, ^4) 



pq 



dt5Gip{ti, is) 



'X-G jq{t2,t^)^pqkl{t^t3) 



(153) 



Because of the equal-time delta function in Eq. (|153() we 
find that 

tr{A"} = 

= j d{ll' . . . nn'){n'\A\22') . . . {nn'\A\ll') 

= / d{ti . . .tn)Ap-^p^{ti,t2) ■ ■ ■ Ap^^p-^{tn,ti) 

^ Jo 



— lim > e^^ Ap-^p^{z) . . . Ap p-^^{2 

'q^0+ 



(154) 



where pfe — (ikjk) are multi-indices and where we defined 

Aijkiititz) = / dt<^Gip{ti,t<^) 

pq Jo 

xGjq{tl,t5)jpqkl{t5t3) (155) 

which in frequency space attains the form 



Aijkl{z) = 4 X! Gip{zi)Gjq{z - Zl)^pqkl{z) (156) 
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From diagonalization of Apq{z) where p = (ij) and q = 
(kl) we then immediately obtain 

Lpp[G, r] = 

Y^iHl + Ap(^)) - Ap(z) + ^Xliz)) (157) 

z,p 

where Xp(z) are the eigenvalues of A{z). Let us finally 
concentrate on the evaluation of B. This expression is 
given by 

Bijklitltltzti) — Gqi{t4, tl) 

PI 

y.Gjp{t2,t'i)-ipMq{hti) (158) 

This expression depends on three relative times which 
makes it awkward to evaluate the logarithm. We there- 
fore follow references^ and replace in frequency space 
li]ki{z) by its zero-frequency hmit "ftjki{0), 

li]ki{hti) ^ ltjki{0)S{t3,t4) (159) 

such that 

Bijkl{tlt2t3t4) — (5(^3^4) Gqi (t4, tl) 

pq 

XGjp{h,t3)-fpklq{0) (160) 

Then, similarly as for the quantity A we have 

tr{S"}= lim ye'^'Bp,pJz)...Bp^pJz) (161) 

where 

B^Jkl{z) ^^Y. Gq^{zi)Gjp{zi + z)-1pklq{Q) (162) 

Now B{z) is readily diagonalized with respect to its two- 
particle indices to give 

Lph[G,f] = 

Y,{\ ln(l - Ap(z)) + \\p{z) + \\liz)) (163) 

where \p{z) are the eigenvalues of B{z). The full func- 
tional L[G', r] is then constructed as 

L[G, f ] = Lpp[G, f] + Lpn[G, f] - ^tr {A^} (164) 

o 

where the last term is easily found by summing the 
squares of the eigenvalues of A and performing a fre- 
quency sum. It finally remains to calculate an explicit 
expression for E[G, F] and to evaluate the last term 
in Ea. p44|l . The self-energy is readily calculated from 
Eqs.lO and to be 

E.,(z) = S,^^ + S,,-c(^) (165) 



where 

S,^^ = \ hm ^ Vo,,pq,Gqp{z) (166) 

z 

and 

S.y,c(2) = ^ 51 yQ,ipqTGrs{z\)Gqt{z2) 
zi,Z2 pqrstu 

xGup{zi +Z2- z)7st«j(0) (167) 

This expression is, of course, considerably simplified 
when we use a diagonal input Green function. This fi- 
nally concludes the discussion on the practical evaluation 
of the functional. 

In summary: evaluation of the functional Y in prac- 
tice therefore essentially involves the diagonalization of 
the one-particle matrix A{z) of Ea. H143|l and the diago- 
nalization of the matrices A{z) and B{z) of Eans. (|155|l 
and l|162(l in a two-particle basis followed by a frequency 
summation. This is, for instance within the DMFT ap- 
proach used by Katsnelson and Lichtenstein22, a numer- 
ically quite feasible procedure. 

IX. CONCLUSIONS 

In this work we studied variational functionals of the 
Green function and the renormalized four-point vertex 
in order to calculate total energies for strongly corre- 
lated systems. The variational functionals were derived 
by Legendre transform techniques starting from an ex- 
pression of the action (or grand potential) defined on the 
Keldysh contour. The structure of the functionals was 
further analyzed by means of diagrammatic techniques. 
We finally gave a detailed discussion of the practical use 
and evaluation of these for different approximate func- 
tionals. Future applications along the lines described are 
intended. 

Finally we comment on further applications of the vari- 
atonal functionals. It was found that within the $ and 
the ^'-formalism could be succesfuUy used to derive ex- 
pressions for response functions within time-dependent 
density- functional theory (TDDFT)^*. This was done by 
inserting approximate Green functions G[v], coming from 
a noninteracting system with a local potential v, into the 
variational functionals. Then the potentials were opti- 
mized by requiring that dY/6v = 0. Due to the one-to- 
one correspondence between the density and the poten- 
tial (as follows from the time-dependent generalization 
of the Hohenberg-Kohn theorem.^) this then implies that 
we are optimizing a time-dependent density functional. 
The optimized potentials are then to be interpreted as 
Kohn-Sham potentials. In this way one obtains a den- 
sity functional for every diagrammatic expression from 
the <&- or ^f-functional. A similar procedure can now be 
carried out for the S-functional. 

A further point of future investigation is concerned 
with finding the variationally most stable functional. It 
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was already mentioned that the Klein and Luttinger- 
Ward (LW) forms of the functional lead to different re- 
sults. The Luttinger-Ward form was found to be more 
stable. This is probably due to the fact that the second 
derivatives of the LW functional are smaller than those 
of the Klein functional. However, it is very well possi- 
ble that one could derive a better functional that would 
make the second derivatives even smaller or make them 
vanish. In that case the errors we make would be only 
to third order in the deviation AG of the input Green to 
the true self-consistent one. This still remains an issue 
for future investigations. Finally we mention that work 
on implementation of the formalism discussed here is in 
progress. 



Then the one-particle Green function in the presence of 
the nonlocal field u is defined as: 



G„(ll') 



5Y 



1 



5Z 



5u{Vl) Z[u] 5u{Vl) 



(A5) 



When evaluated at u = the Green function reduces to 
the familiar one 



G„=o(ll') = -*(Tc[^J,(l)^(l')]) 



(A6) 



Let us note that one should be careful with dealing 
with time-nonlocal potentials. It would, for instance, be 
tempting to think that G„ would be given by the expres- 
sion 
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APPENDIX A: A GENERATING FUNCTIONAL 
FOR THE GREEN FUNCTION 

In order to obtain the Green functions as variational 
derivatives we define^! an evolution operator in terms of 
a time- and space nonlocal potential w(12) 



C/M(to - i/3,to) = 7bexp(-i j dtH{t)- 
i [ dl [ d2V;t(xi)u(12)V;(x2)) (Al) 



where we used the compact notation 1 = (xi,ti). Since 
we have now two times in the exponent this expression 
only has meaning if we define how the time-ordering is 
specified if we expand this expression. It is defined as 
follows: 

U[u]{tQ - t(3,to) = U[u = 0](to - il3, to)+ 

{-^y 



E 

11=1 



d{n' . . .nn') u{l'l) . . .u{n'n) 



x(Tc[V31^(l')V^ff(l)...V'l^(n')^HH]) (A2) 

where the expectation values under the integral sign are 
averages (as in Eq.©) in the absence of the nonlocal 
field u. This definition agrees in the limit of a time- 
local potential, i.e. u(12) = u{jiiti,JC2ti)6{t^ ,t2), with 
an expression that can be derived directly from the time- 
dependent Schrodinger equation. We now define the 
functionals 

Z[u] ^ Tr{U[u]ito-il3,to)} (A3) 
iY[u] = ^ In Z[u] (A4) 



G„(ll') 
Tr 



[u[u]{to-tf3,to)Tc[i^l,{l')i^H{l)]} 
TT{U[u]ito-if3,to)} 



(A7) 



where the Heisenberg operators in the presence of u 
would be given by Oh — U[u]{to,t)OU[u]{t,to)- How- 
ever, this expression is not valid if u is nonlocal in 
time. For instance, when expanding Ea. (|A5p and (|A7|) 
in powers of u using Eq. (|A2|I one immediately sees that 
certain time-orderings of the field operators in Ea. (|A5p 
are absent in Ea. ljA7p . Similarly the evolution operator 
U[u](t^t') does not satisfy a simple equation of motion. 
However, we can still derive the equations of motion for 
Gu on the basis of the hierarchy equations of the n-body 
Green functions in the absence of the nonlocal field u, as 
we will show below. 

More generally we can now define n-body Green func- 
tions Gn,u from a repeated differentiation of Z[u\, i.e. 



Z[u] SuiVl)...Su{n'n) 



iG„,„(l...n, l'...n') 



(A8) 



where e„ = (_i)n("+i)/2_ -phe prefac- 

tor e„ results from reordering the opera- 
tor product {Tc['4>''{l')i){l)...tp^n')i>{n)]) to 
{Tc[ip{l) ■■ ■'ilj{ri) if}'' {!').. .^p''{n')]) as is easily veri- 
fied by induction. One can readily check that for u = 
Eg. ljASp agrees with our previous definition of the 
n-body Green function of Ea. H19|l . From Ea. (|A8|) we 
further immediately obtain that 



1 



(5m(32) 
SZ SZ 



(- 



1 



SZ 



Z[u] 5u{Al)' 
1 5'^Z 



Z2 5u{U) (5u(32) Z 5u{Al)6u{?,2) 
= G„(14)G„(23)-G2^„(1234) (A9) 

This is Eq.|2Zjl used in section lTTll As a next step we de- 
rive the hierarchy equations for the Green functions Gn,u- 
From Ea. (|A8p we see that we can expand the functional 
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Z[u\ as a Taylor series expansion in u as 

ZM = Z[0]^^ / d{ll'...nn') 

xG„(l . . . 71, 1' . . . TiXl'l) . . . u{n'n) (AlO) 

where the term with n = is just defined to be one. 
The Green functions Gn are the Green functions in the 
absence of the field u. The one-body and n-body Green 
functions can therefore be expressed in terms on the field- 
free Green functions using Eas. (|A5p . HA8|I and IjAlOp . 
One obtains for Gu and 6*2,11 the equations 



Z[u] 

m 



Z[u] 

m 



G„(11') = G(11') 

xG„(l . . . 71, 1' . . . 7i')u(2'2) . . . u{n'n) (All) 
G2,„(12;1'2') = G(12;1'2') 



E 



{n - 2)! 



n— 3 



d(33' . . . nn') 



X G„(l . . . 71, 1' . . . 7i')u(3'3) . . . u(7i'7i) (A12) 

From these expressions we see that G„ and G2,u inherit 
the Kubo-Martin-Schwinger boundary conditions from 
the Gn- If we act with the operator idt^ — h{l) on both 
sides of Ea. (|All|l and use the Martin-Schwinger hierar- 
chy equations Ea. H21l) for the Green functions G„ in the 
absence of the u-field together with Ea. (|A12|l we obtain, 
after slightly tedious but straightforward manipulations, 
the equation of motion for Gu 

{idt,~h{i))Gu{n')^5{n') 

d2u(12)G„(2l') 



-i j dxi;(xi,x)G2,„(l,xii,xt+,l') (A13) 

By functional differentiation with respect to u we can 
generate equations of motion for the higher-order Green 
functions. To see this we first multiply <jA13|) by Z[u] to 
obtain 

{idt,-h{i)){-i){\i') = 5{n')z[u] 

+ j dlu{ll){-i){ll') 

-i j dx7;(xi,x)(-i)2(l,xti,xij^, 1') (A14) 

where we introduced the simplified notation 

(1 . . . n; 1' . . . n') = Z[u] G„,„(l . . . ?i, l' . . . n') (A15) 

We use the convention that the primed variables are al- 
ways associated with creation operators and that the 



unprimed variables are always associated with the an- 
nihilation operators. Taking the functional derivative of 
Ea. (jA14p with respect to u(2'2) then gives 

{idt, - ;i(l))(-i)2(ll'2'2) = 5(ll')(-*)(2'2) 
+5(12')H)(21') 

+ j dlu{ll){-if{ll'2'2) 

-i J dxi;(xi,x)(-i)^(l,xii,xi+,l'2'2) (A16) 

Reordering the indices and dividing by Z[u\ then gives 

(ia*, -/7(1))G2,„(121'2') = 

-(5(11')G„(22') + (5(12')G„(21') 



<ilu(ll)G2,„(121'2') 
-i y"dxt;(xi,x)G3,n(12,xti,x<+,l'2') (A17) 

By continued differentiation we obtain the general hier- 
archy equations for Gn.u 

-/i(l))G„,„(l...n,l'...n') 

n 

= ^ <5(lj')(-l)"-^G„_i,„(2 . . . n, 1' . . . / - 1, j' -f 1 . . 



+ J dlu(ll)G„,,(12...71,l'...7l') 

-i J dxv{xi,x.)Gn+i,ui^ ■ ■ ■n,xti,xtf ,1' . . .n') (A: 

These equations are readily checked by induction if we 
multiply them by Z[u\ and take the functional derivative 
with respect to 71(71' -I- l,n -I- 1). We have therefore es- 
tablished that Green functions Gn,u satisfy an obvious 
generalisation of the hierarchy equations. The relation 
Ea. (|A18|l is the main result of this Appendix and will 
be essential to the derivation in the next section. Note 
further that Ea. HA18|l can be used to derive a Wick's 
theorem in the presence of the nonlocal field u. If we 
put w = we find that the noninteracting 71-body Green 
functions G„,„ satisfy Ea. (jA18|l if they are written as 
determinants in terms of Gu- 

APPENDIX B: THE EQUATION OF MOTION OF 

Gu,v 

The main goal in this Appendix is to derive the equa- 
tions of motion for the Green function in the presence of 
the a nonlocal one-body potential u(12) and a nonlocal 
two-body potential 1/(1234). As in Appendix IXIthe main 
difficulty is caused by the fact that u and V are nonlocal 
in time. Our final result can be obtained with help of 
Ea. (|Il8)l . Let Z[u, V] be given by 



Z\u,V] - 



V]^Tr{u[u,V]{to'tP,to)} 



(Bl) 
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where 



U[u,V]{to - il3,to) = TceM-i / dtHo{t) 



i I dl J d2V't(xi)u(12)7/'(x2) 
d(1234)F(1234) 



XV'^(X1)V'^(X2)V'(X3)V'(X4)) 



(B2) 



Due to the multiple time-integrals this expression has 
only meaning when we define how the time-ordering is 
specified when we expand this expression. We define 

U[u, V] {to - ^/3, to) = Lf[0, 0] {to - if3, to)+ 



E 



n 



!to!4 



ni 



divi ■ ..yn)d{Xi . ..Xrn) 



n,m— 1 

XMi . . . UfiVi . . . XmTc[Vl ■ ■ ■ VnXi . . . X„i] (B3) 

where for the coordinates we introduced the short nota- 
tion 



K 
X, 



(z'J) (B4) 

((2i- l)',(2i)',2i- l,2i) (B5) 
and we further defined 

U], = u{y,) (B6) 

^P^Hii')^H(i) (B7) 



V{X, 



(B8) 



X, = ^^((2i - l)')V3Ij((2z)')^H(2i - l)^H{2i){m) 

where the Heisenberg representation of the operators 
is defined with respect to Ho- We then define the n- 
particle Green function Gn,u,v in the presence of the 
time-nonlocal fields u and V as 

enG„,„,y(l . . .n; 1' . . .n') = 

1 ^"^[^,^] , . 

Z[u,V]5u{Vl)...u{n'n) ^ ' 

where e„ = (-l)"("+i)/2. For = this definition 
agrees with the definition (|A8() in Appendix El in the 
absence of two-particle interactions. Now from Ea.(|B3|l 
one can readily derive that 

5^Z[u,V] 



5Vi... 5Vk 



S^''Z\u,V] 



4fe 6u{l'l)...u{{2ky,2k) ^^^^^ 

With this equation we find that we can express the Green 
functions equivalently as 

G2„,„.v(1...2n; l' . . . (2n)') 

4"(-i)" S''Z[u,V] 



Z[u, V] SVi... SVn 
G2„+i,„,y(l,l,...2n;r,l'...(2n)') 

" Z[u,V] 6uiVi...Vn 



(B12) 



(B13) 



As a particular case we have 
G„,y(ll') = 



1 SZ[u,V] 
Z[u,V] Su{l'l) 



(B14) 



These equations are two basic starting equations Eq.( 
I41|l and 1)42(1 of section Hvl It remains to show that they 
are related by an equation of motion. From Eas. l|B12|) 
and HB13|I we see that Z[u, V] has the following Taylor 
expansion around V = 0. 



Z[u, V] _ 
Z[u,0] f-n!4 



d{X,...Xn) 



G„,y(il') = G„(1,1') 



xG2„,„(l . . . 2n; 1' . . . 2n')Vi . . . K (B16) 

where we denote Gn,u = Gn,u,v=o- Similarly for Gu,v 
we have from Ea. ((B13|l the Taylor series expansion 

Z[u,V] 
Z[u, 0] 

+ E;^/^(^i---^«) 

n=l ■ 

xG2„+i,„(ll ...2n; I'l' . . . 2n')Vi . . . K (B17) 
and for the two-particle Green function from Ea. l)B12|l 
Z[u,V]^ 



Z[u, 0] 



G2,n,y(12,l'2') =G2,„(12,1'2') 



- ^Ei^I-T^/'^(^i--^-) 

n=2 ^ ' '' 

xG2„,„(1212 ...2{n- 1); r2'l'2' . . .2{n - 1)') 

xFi...K-i (B18) 

To obtain an equation of motion of Guy we can use the 
hierarchy equations for Gn,u of Eq. for ti; = 

-Ml))Gn,n(l...n,l'...n') 
xG„_i,„(2...n,l'...(j-l)',0' + !)'...«') 



+ / dlM(ll)G„,„(12...n,l'...n') 



(B19) 



If we act with idtj — h{l) on both sides of Ea. l|B17|l and 
subsequently use Eas. ((B19|l . l(B16|l and ((B18|l we obtain 



{idt,-h{i))Gu,v{n') = sin') 

d2u{12)Gu,v{2l') 



--J d(234)F(1234)G2,„,y(432r) (B20) 

This is the equation of motion for the Green function 
used in section llVl Again by differentiating this equation 
with respect to u we obtain the hierarchy equations for 
the higher order Green functions Gn,u,v- 
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APPENDIX C: FEYNMAN RULES FOR THE 
TWO-PARTICLE GREEN FUNCTION 

In this section we give a brief summary of the Feyn- 
man rules for the two-particle Green function G2 within 
the Hugenholtz diagram techniqu6SS*2SiSL2&. The general 
structure of the two-particle Green function is as given in 
Fig.©. The Green function 6*2(1234) is written with the 
points (1234) positioned clockwise on four corners of the 
diagram where corners 1 and 2 are connected to outgoing 
lines and corners 3 and 4 are connected to ingoing lines. 
If one expands the evolution operators in the definition 
of G2 in powers of the interaction V one finds for the 
diagrams the following rules 

1. Every Green function line (contraction according 
to Wick's theorem) gives a factor iG . 

2. Every vertex gives a factor —iV . 

3. Every closed loop of Green function lines gives a 
minus sign, i.e. we have a factor (—1)' where I is 
the number of closed loops. To find the number of 
loops one must replace the Hugenholtz vertex by 



the first term on the right hand side of Fig|^ (with 
the same labeling) and count the number of loops 
that appear in this way. 

4. A line starting at 3 and ending at 1 gives a minus 
sign, i.e. we have a factor (—1)^1^ where L13 = 1 
when 1 and 3 are connected and zero otherwise. 
To determine the connectivity it is necessary that 
one again first replaces the Hugenholtz vertex by 
the first term on the right hand side of Fig|21 (see 
also21). 

5. Two Green function lines (so-called equivalent 
lines) starting from a given vertex and ending both 
on the same vertex give a factor i, i.e. we have 
a factor 2~p where p is the number of equivalent 
lines. 

6. There is a factor (— i)^ from the definition of G2. 

From these rules we find that the overall prefactor of a 
G'2-diagram with n vertices is given by j"(— 1)'+^"2~p. 
For example, diagrams (a) — (/) in Figl^have prefactors 
1. — 1, i, — 1 and —1 respectively. 
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